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Abstract 

We study the discrete massless Gaussian Free Field on Z d , d > 2, in the presence of a disor- 
dered square-well potential supported on a finite strip around zero. The disorder is introduced by 
reward/penalty interaction coefficients, which are given by i.i.d. random variables. 

In the previous note [J], we proved under minimal assumptions on the law of the environment, 
that the quenched free energy associated to this model exists in K + , is deterministic, and strictly 
smaller than the annealed free energy whenever the latter is strictly positive. 

Here we consider Bernoulli reward/penalty coefficients b ■ e x + h with e x ~ Bernoulli 1/2 +1) 
for all x € Z d , and b > 0, h € M. We prove that in the plane (6, h), the quenched critical line 
(separating the phases of positive and zero free energy) lies strictly below the line h = 0, showing 
in particular that there exists a non trivial region where the field is localized though repulsed on 
average by the environment. 
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1 The model 

We study the discrete Gaussian Free Field with a disordered square-well potential. For A a finite 
subset of Z rf , denoted by A <s Z d , let cp = ((p x )xeA represent the heights over sites of A. The values of 
(p x can also be seen as continuous unbounded (spin) variables, we will refer to ip as "the interface" or 
"the field". 

Let Q = M z be the set of configurations. The finite volume Gibbs measure in A for the discrete 
Gaussian Free Field with disordered square- well potential, and boundary conditions, is the probability 
measure on Q defined by : 

MA°( d V) = 7^0 6X P ( -^A(f) + P ' ex + /l ) 1 [^6[-a,a]] J n dip x Yl $o(d<p y ). (1) 

A \ xeA / xeA yeA c 

where a, (3,b > 0, h E M and T~L\((p) is given by 

U ^ = T d E fa-Vv) 2 ' ( 2 ) 

{x,y}nA^0 



where x ~ y denotes an edge of the graph 7L d and Ir^i denotes the indicator function of A. An 
environment is denoted as e := (e x ) xe \. We consider here e given by i.i.d. random variables 

e x ~ Bernoulli 1/ / 2 ( — 1, +!)• 

The parameter b is usually called the "intensity of the disorder", while h is its average. The disordered 
potential attracts or repulses the field at heights belonging to [— o, a]. Z^f is the partition function, 
i.e. it normalizes /j^f so it is a probability measure. The superscipt reminds the boundary condition, 
it is added to the notation compared to [4J because it will be useful below. We stress that our model 
contains two levels of randomness. The first one is e which we refer to as "the environment". The 
second one is the actual interface model whose low depends on the realization of e. 

The inverse temperature parameter (5 enters only in a trivial way. Indeed, if we replace the field 
(<Px)xeA b y {VP<f>x)xeA, a by ^a, and (b ■ e x + h) xeA by (J3(b ■ e x + h)) x( z A we have transformed the 
model to temperature parameter (3=1. In the sequel, we will therefore work with f3 = 1. 

The dimensions 1 and 2 are physically relevant as interface models. In this paper we focus on 
d > 2 since 1-dimensional models have been well-studied in the last decade (see |lj for a historical 
introduction). 

The questions we are addressing in this framework are the usual ones concerning statistical me- 
chanics models in random environment : Is the quenched free energy non-random ? Does it differ from 
the annealed one ? Can we give a physical meaning to the strict positivity (resp. vanishing) of the 
free energy ? What can be said concerning the quenched and annealed critical lines (surfaces) in the 
space of the relevant parameters of the system ? 

In the previous note [3J, we proved under minimal assumptions on the law of the environment, that 
the quenched free energy associated to this model exists in M + , is deterministic, and strictly smaller 
than the annealed free energy whenever the latter is strictly positive. 

Here we investigate the phase diagram of the model : in the plane (b, h), we prove that the quenched 
critical line (separating the phases of positive and zero free energy) lies strictly below the line h = 0. 
Thus there exists a non trivial region where the field is localized though repulsed on average by the 
environment. 



2 Results 

We define the quenched (resp. annealed) free energy per site in A <s Z rf by : 

/ 7 e,o\ /F7 e '°\ 
#(e) = IM- 1 log ( J , /1(e) = lAr 1 log ^ , (3) 

where denotes the partition function of the model with no potential, e x = (i.e. of the Gaussian 
free field). In the case when A = A n = {0, n — l} d we will use short forms f% (e) and /^(e). By the 
Jensen inequality, we have / q (e) < / a (e). Moreover, it is not difficult to see that the annealed model 
corresponds to the model with constant (we will also say homogenous) pinning with the strength 

1(e) := log(E(e b - e * +?i )). (4) 

for all x G A. In other words EZ^'° = . 

In [U Fact 2.3] we proved that for any environment e such that the annealed model exists, i.e. 
K(e b ' ex+h ) < oo, both the quenched and annealed free energies are non-negative. This motivates the 



2 



following notions. We introduce the quenched (resp. annealed) critical lines, which are delimiting the 
region where / q (e) = (resp. f a (e) = 0) from the region / q (e) > (resp. / a (e) > 0). 



h%(b) := sup{/i € R : / q (e) = 0} and h*(b) := sup{/i G R : / a (e) = 0} 

We are interested in describing the behavior of these quantities in the phase diagram described by 
the plane (b,h). Knowing the behavior of the homogenous model for positive pinning [3], we easily 
deduce that the annealed critical line is given by the equation 1(e) = 0. 



Note that / q (e) < / a (e) implies that h^(b) > h*(b). In Theorem 2.1 we show that the quenched 
critical line lies strictly below the axis h = in the neighborhood of b = for all d > 2. Our result 
shows in particular that there exists a non trivial region where h < 0, b > and / q (e) > 0, i.e. where 
the field is localized though it is repulsed on average by the environment. 
Note that we don't have any estimate on the behavior of /i q (6) — h*(b). 

Theorem 2.1. Let e ~ ® X £ Z dBernoullii/2(— 1, +1)- Then, 

For d>2, the quenched critical line is located in the quadrant {(b, h) : b > 0, h < 0} . 

More precisely, there exists some C, C > depending on d, a only and e G (0, 1) such that for any 
environment e which fulfills b + h > ; — e < —b + h < and 



h > -C ■ (-b + h) A ford>3 



we have / q (e) > 0. 



Remark 2.2. 1. A sketch of these bounds in the plane (b,h) can be seen on Figure [7J Moreover, 
the bound for d > 3 can be rewritten as h > —C"(d, a) ■ b 2 . 

2. Jensen's inequality gives us an upper bound on C, C . Indeed, as / a (e) > / q (e), if / a (e) = 
then / q (e) = 0. In particular, we must have —C < j^h c \ b= Q < 0. Our result gives thus an 
upper-bound on the behavior of the quenched critical line near 6 = 0. 



2.1 Related results 

The same type of result has been proven for (Id) polymer models in great generality in [I], where 
Alexander and Sidoravicius consider a polymer, with monomer locations modeled by the trajectory of 
a Markov chain (Xj)j e g, in the presence of a potential (usually called a "defect line") that interacts 
with the polymer when it visits 0. Formally, the model is given by weighting the realization of the 
chain with the Boltzmann term 

exp (^Y^(u + Vi)t [Xi=0 ^j . 

with (Vi)i£z an i-i-d. sequence of 0-mean random variables. They studied the localization transition 
in this model. We say that the polymer is pinned, if a positive fraction of monomers is at 0. In the 
plane ((3,u) critical lines are defined as above: for f3 fixed, u^((5) (resp. u a (/3)) is the value of u above 
which the polymer is pinned with probability 1 (for the quenched (resp. annealed) measure). They 
showed that the quenched free energy and critical point are non-random, calculated the critical point 
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Figure 1: (Color online) Phase diagram of the model. The red curve is the annealed critical line; the 
blue one is our bound on the quenched critical line. 



for a deterministic interaction (i.e. Vi = 0) and proved that the critical point in the quenched case is 
strictly smaller. 

When the underlying chain is a symmetric simple random walk on Z, the deterministic critical 
point is 0, so having the quenched critical point u c {(3) strictly negative means that, even when the 
disorder is repulsive on average, the chain is pinned. This result was obtained by Galluccio and Graber 
in [7] for a periodic potential, which is frequently used in the physics literature as a "toy model" for 
random environments. 

A much shorter proof with explicit bounds can be found in [8], in less generality, but [6] contains 
a revisited proof with explicit estimates and weakening the assumptions on the underlying model. 

Note that for polymers, or discrete height interfaces, one need a coarse graining procedure to achieve 
the proof. In our case, as we will see in the next section, we can shift the continuous interface where 
the environment is unfavorable, and this has a small cost in dimension d > 3. The procedure is a bit 
more complicated in dimension 2 and we have to localize the field before by introducing a small mass. 

2.2 Proof of Theorem O 
2.2.1 Cased>3 

We assume —b + h<0<b + h i.e. the environment is repulsive if e x = — 1 while it is attractive if 
e x = +1. 

The idea is to tilt the measure such that the field ip is shifted up of an amount s on the sites x 
for which b ■ e x + h < 0. In this way the shift of the field follows the environment. For some technical 
reasons, we need to work with the measure with boundary condition a, so perform two changes of 
measure (first one changing boundary condition and the second one to following the environment). Let 
s > (to be fixed later). 

(, o,o , o,o / 
~ A Z n ft £ ^ ■ e * + Veha,a]] 

where Wx) x& x n under ^An.e.s is distributed as (ip x + st[( b . ex +h)<0]) xeAn under Ma*- More formally, 
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introducing T e , s : {((p x ) xeAn ) ^ {fx + sl l(b-e x +h)<o]) xeAn , we define as n£ o T e J. 

Using Jensen's inequality, we get 



fUe) = n- d log 



0,a 



j 0,0' 



• in ' i 0,a 



^ — d 0,a 
^ n »A n ,e,s 



i 0,a i 0,0 

E (b ■ e x + ^)l[^ e[ _ a , a]] + log + log -% 

. * v ' * v '. 

V (1) (2) / 



As Z 



0,0 



, , c s — Z^® (which follows by change of variables in the Gaussian integral), the first term can 
be written as 

W = ~AJ E ~Vyf ~ (<Px -<Pyf 



{x,y}nA„^0 

x~y 



where tp x := tp x + sl[ b . ex+h<0 y Hence, using the definition of Marcs' 

/ 



n 0,a 

" Ad ^ 



E - ^v) 2 ~ (V>x ~ <Py) 



{x,y}nA n ^0 

\ x~y 



) 



s 2 n d 
Ad 



E ( 1 [b-e x +h<0] ~ ^[b-e y +h<0]Y 



{x,y}nA n ^0 
x~y 



The second term contains only boundary contribution of order n d 1 . Indeed, 



Z 



0,cT 



(2) = 2a £> x - « 2 l aA «l + lo S dfe ^ 2a E V* ~ fl2 l aA «l " Cn ' 



d-i 



xedAn 



xedAn 



Hence, 



n 



"V A ': e , s ((2))>2a-n- rf £ M A '>*) ~ ^n" 1 > s E W*<o] - Cn~ l > -Cn' 1 



xedA n 



xedA n 



We get 

fl(e) > n- d Y,(b-e x + h)^ x G[-a,a]) 



xeA Tl 



s 2 n d 
Ad 



E (Mb-ex+hKO] ~ l[6-e„+/i<0]) 2 ~ C n 1 



{x,y}nA„^0 



Now we use the fact that the marginal laws of all ip x , x G A n under are Gaussian variables centered 
at a, i.e. ip x ~ Af(a, ) where cr£ = Var A ' a (<^) < Var^fax) < c(d) < oo for d > 3. Therefore, 



Ma^P* g [-a,a])-^(f x + se [-a, a}) = n u £((p x G [-2a, 0]) - f^(<p x G [-2a-s,-s]) 



0,a 



0,0/ 



0,0, 



-2a 



2a— s 



s as n — ^ oo, 



(5) 
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for c(d) >> s. In particular we will use that : 

for some C\(d, a) > 0. 

#(e) > n~ d £ (6 • + (^(^ G [-a, a]) - d(d, a)st [b . ex+h<0] } 



xeA n 



s 2 n d 



^ (l^-es+ZKO] - l[6-e H +/i<0]) 2 ~ Cn 1 



Ad 

{x,y}nA n7 ^0 
xr^y 

Observe that [j^ a (tp x 6 [—a, a]) = /J^°(ip x G [— 2a, 0]) > Lk£(ip x £ [— 2a, 0]) > C2{d,a) for some 
C 2 (d,a)>0. 

By taking the expectation with respect to the environment, using the bounded convergence theorem 
and the fact that /1(e) = E(/ q (e)) (cf.0| Theorem 2.1]) we get : 

/»=JimE/^(e) > h C 2 (d,a)- SC ^ a ^(-b + h)-^ (6) 



n— >oo 



We may optimize over s as the left hand side does not depend on it. Doing this one checks that 
/ q (e) > as soon as 

C2{d,a) 

This gives the implicit equation in terms of the variance b 2 of b • e x + h : 



h>b- — + 1 J^-- = -Kb 2 + 0(b 3 ) 
2K 2 \ K 2 K K ' 

The annealed critical curve as well as this bound are drawn on Figure [T] We recall that ^ is valid 
under assumption that s is small. The maximum of ^ is realized at s max = — 4Ci ■(—& + /&), thus it 
is enough to assume that (— b + h) is small. □ 



2.2.2 Cased = 2 

In the case d = 2, the variance of the Gaussian free field diverges with the size of the box, so we cannot 
use the previous estimates. To circumvent this problem we introduce the so-called massive free field. 
Let m > 0, 

vY n ,m( dl P) = exp ( ~ m2 Yl &x ~ C) 2 ) 11 d<p x Yl Soi&Px), ( 7 ) 

A n ,m \ xeA n / xeA n x£dA n 

where T-L\{<p) is defined in ([2]). Known facts about this model can be found in [5j Section 3.3]. In 
particular, the random walk representation for the massless GFF j3j (1.3)] is still true, but for a 

2 

random walk Yt that is killed with rate £(m) = 1 ^ Tt a , namely at each time £, if the walk has not 
already been killed, it is killed with probability £(m), where the killing is independent of the walk. We 
write its law P x when it starts at x. 
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Lemma 2.3. Let d = 2. Then, 

1. There exists some C\ > such that for n large enough, m > small enough and all x S A n , 

/xSf )m (^)<Ci|Iog(m)|. 

2. There exists some C2 > suc/i i/iaf for n large enough and m > small enough, we have 

Z°'° 

n~ 2 log A q'^ > -C 2 m 2 | log(m)|. 

Jin 

Proof. These bounds are rather standard. We give here the main steps of the proofs with some 
references. For the first claim, we use the random walk representation [5] to write 

00 00 

where ta„ is the first exit time of A n and H is the killing time of the random walk Yt. Hence, 

00 

l*Z,J&) ^ »Z,m(<Pl) < E(l " tdn)fP»{X t = 0) (8) 

where Xg_ is a simple random walk (without killing). The projections of X(_ onto the two coordinate 
axis are two independent 1— dimensional random walks Xj and Xf, then by Stirling formula, 

Po(x 2i = 0) = (p (x l 2e = o)) 2 = (( 2 /V 2 *V = + o(1)) as e ^ 00 

The asymptotics of Q for small m gives the desired upper-bound. 

To prove the second claim we use the representation of the partition function described in [5J p. 542] (it 
applies to the massive GFF with an obvious modification). We denote by P the coupling of a random 
walk X n and a killed random walk Y n such that Y n = X n up to its killing time H, 



Z°'° 

A n ,m 



V x£\ n e=i 



Z°'° 

l A «r llo g^o4^ = l A nl~" U > I > ^ = x ^A n > 21) - P x (Y 2i = x,T An A N > 2£) 



1 A 1 



00 



11 

^-F o (I a = 0,TA B >2£,K<2£) 



2^2^ 
1" U 



^ 2 E/ofe = 0)(l-(l-rf) (9) 

£=1 

Using the same estimate as in Q, the asymptotics of ^ for small m gives the desired upper-bound. □ 

The idea is to tilt the measure, as in the proof for d > 3, first to work with the massive measure, and 
second to follow the environment such that the field (p is shifted up of an amount s on the sites x for 
which b ■ e x + h < 0. For some technical reason, we need to work with the measure with boundary 
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condition a, so we perform three changes of measure (first one for changing boundary condition, a 
second one for adding mass, and a third one for following the environment). 
Let s > and m > to be fixed later. 

./A n ( e ) = n ~ 2 iog^A'! ( ex P ^2(b-e x + h)t [lj)xe[ _ a>a]] J 

= n- 2 log/x^ m e)S I exp ( J2 ( b ■ e * + /l ) 1 [0,e[-a,a]] 
V \xeA n 



a 0,0 



+ log T-oZ +log A o,c 



d/J°' C d/z°' C m 

' v - +log- 



du°' C 

A n ,m "A n , 771,6,5, 



where (^x) :ce A n under /*Slm,e,« is distributed as + «l[(6. eit +h)<o]) a . 6AB under t*j£, m 5 more formally, 

introducing r e>s : ((^) xeA J (<p x + sl [(&-e,+fc)<o]) a . eAn . we define /4tm,e,s as /z°'f m o T^ 1 . Using 
Jensen's inequality we get 

/ \ 

o,o j.o.c 



fl(e) > „-VK 



j o,c 



1 U,U j I 

£ ( & • e * + ^)l[^6[-,a]] +'°gT^ + log TlSr- + bg ^ ■ 

d MA n d^At.m d^A^.m.e.s 



xGA 



o,C 



V 



(1) 



(2) 



(3) 



/ 



As in the proof for d > 3, we have 



By Lemma 



2.3 



we have 



r0,f 



7 o,c 7 o,o 7 o,( 



A n ,m A 



tr^j — ox — —Cn — C2TI m \ logm|, and then 

Z\ z 



/ z o,C \ 

(2) = log -^p \ +m 2 vl> ~ Cn ~ C 2 n 2 m 2 \ logm| + m 2 ^ ip% 

V ^A n / a;GA n x£A n 



hence, 



n 



Finally, noticing that Z ^ mes z 
we can compute the third term. 



V A f ,™,e, s ((2)) > "Cn" 1 - C 2 m 2 | bgm|. 

(just perform a change of variables in the Gaussian integral), 



z o,C 

A„,m 



(3) = -J £ (Vx ~ <Py) 2 ~ (fix ~ fiy) 2 ~ m 2 ^ (<p x - sf - (<f x - S)' 



{x,y}nA n ^0 
x~y 



xeA, ; 



where <p x := cp x + sI^.^+zko]- Now we will use the fact that the marginal laws of all <p x , x S A n 
under are Gaussian variables, i.e. ip x ~ Af '(/i^, c^ 2 ) where /i^ ~ C except for x close to the 

boundary of the box. Indeed, by the random walk representation of the mean, there is C > such 
that |/x°'f jm (^) -C| < C{l + m 2 )-^- dA ^. Moreover, (o x n f = Vai°£ m ((p x ) < Ci|logm|. Using the 
definition of /x^ mes , an d computing the terms as in the proof for d > 3, 



*-V£, w ((3)) * 



S 

8n 2 



( 1 [b-e x +h<0] - 1 [b-e y - 



■h<0], 



{x,y}r\A n ^0 
x^jy 



^ 2 m 2 s 2 
n 2 



C 



£ 1 [b-e.+/ l <0] + — • 



xeA„ 
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We get, for n large enough and m small enough 



ffcje) > n 2 ^ (b ■ e x + h)fx°^ m (0 x G [-a, a]) - ^ ^ {^[b-e x +h<o\ ~ 3 b-e y nj 

x~y 

^ l[5. ex+ ft,<0] - C'm 2 \ logm| - Cn" 1 , 



^2 



m 2 s 2 



n 2 



x<=A n 

,2\ 



for some C and C" > 0. Note that for 0(n ) sites x, we have 

MSlmCv* G h a > a D - MSi, m (V« + s G I _a ' a D ~ $ C,fc( a ) • s as n ^ oo, (10) 

for s << a < C = Ci| logm|, and b = Ci \ logm|. Above stands for the p.d.f. of the above Gaussian 
distribution with mean £ and variance b 2 . In particular, for a positive fraction of x (close to 1) and m 
sufficiently small, we have the upper bound : 

^miVx G [~a, a]) - H^nfa + s G [-a, a]) > • s, 

for some Ci(a) > 0. Now we can compute : 

/^ n (e) > n- 2 ^{b-e x + h)(ji Q £ t Jtp x G [-a, a]) - i J " i gl [^+fe<o]) 

S 2 v a . .2 Til 2 S 2 v a 

t 1 |*-e«+h<0] - 1 [6-e !/ +fe<0]J ~ 2^ 1 [^+^<0] 

{x,j/}nA n ^0 xeA n 
— C"m 2 | logm| — Cn -1 

Observe that /x A '^ m (^x G [—a, a]) > 2a • $^ b2 (— a) = j^gm! uniformly in x G A n . Let us take 
the expectation with respect to the environment, and use the bounded convergence theorem and jH 
Theorem 2.1], we get : 

x , ir~* fO / \ , Ci(a) Ci(a)(— 6 + s 2 m 2 s 2 . 9|1 , . . 
f q (e) = lim E/? e > fc - u ' -s- u , — - Cm 2 log m. 11 



logm| 2|logm 2 16 



Our aim now is to show that the right hand side can be positive even when h is negative. In the above 
expression s, m are free parameters which we may vary. However, we have to remember that both s 
and m need to be small enough, which makes standard optimization analysis cumbersome. We are 
going to show that there exists C > and e > such that for any b, h such that (—6 -+- h) G (— e, 0) 
and 

(-b + h) 2 



h:=C 



log(-(-6 + fc)) : 



there exist small s and m such that the r.h.s. of (11 ) is positive. Notice that the result will imply that 
for any h > C lo ^_^_ b+h ^ the free energy is positive. Let us choose the value of s which maximizes 
(11 ) for fixed m, i.e. 

Ci(q)(-ft + fe) ( , 

(m 2 + 1/4)| log(m)r 1 ' 



9 



and for m let us take 

m 2 = -k/{\ogkf, where k := -hCi{a)/C. (13) 

One can verify that with the above choice of parameters both s and m are as small as we want. Let 
us first put (12) into the r.h.s. of (11) and obtain 

|iogm| (m z + l/4)(logm) z 
For k and consequently m small enough we have 

r e > h- r + — ^ Cm logm. 

|logm| 2(logmj z 

Further let us multiply both sides by (logm) 2 and insert (13). 
/ q (e)(logm) 2 > -/iCi(a)logm+ Cl(a)2( ~ 6 + /t)2 +C / m 2 (logm) 3 

- Qogfc - 31og(| logfcl)) + Cl(Q)2( - 6 + /t)2 - 0*&* k ~ 31 ° g(l l0g " ))3 



2 

Ci(a) 



8(log/c) J 



, , / % Ci(a) 2 (-b + h) 2 . Ci(a) „ 

-/i—^- log(-/iCi (a)/C) + — ^+/i^^+o(/i) as 

2 2 o 

From the last claim it is straightforward to conclude existence of C (sufficiently small) and e with the 
properties described above. □ 
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